Assessment study of ionosphere correction model using single- and multi-shell algorithms approach over sub-Saharan African region

It is a known fact that ionosphere is the largest and the least predictable among the sources of error limiting the reliability and accuracy of Global Navigation Satellite Systems (GNSS) and its regional augmentation systems like Satellite Based Augmentation System (SBAS) in a safety-of-life application. The situation becomes worse in the Equatorial Ionization Anomaly (EIA) region, where the daytime ionization distribution is modified by the fountain effect that develops a crest of electron density at around ±15° to ±20° of the magnetic equator and a trough at the magnetic equator during the local noon hours. Related to this phenomenon is the appearance of ionosphere irregularities and plasma bubbles after local sunset. These may degrade further the quality of service obtained from the GNSS/SBAS system of the said periods. Considering the present operational augmentation systems, the accuracy and integrity of the ionosphere corrections estimate decreases as the level of disturbances increases. In order to provide a correct ionosphere correction to the user of GNSS operating in African EIA region and meet the integrity requirements, a certified ionosphere correction model that accurately characterizes EIA gradient with the full capacity to over-bound the residual error will be needed. An irregularities detector and a decorrelation adaptor are essential in an algorithm usable for African sub-Saharan SBAS operation. The algorithm should be able to cater to the equatorial plasma vertical drifts, diurnal and seasonal variability of the ionosphere electron density and also should take into account the large spatial and temporal gradients in the region. This study presents the assessment of the ionosphere threat model with single and multi-layer algorithm, using modified planar fit and Kriging approaches.     

On the time of minimum ionization in the F2 region

Using measurements of the critical frequency of F2 region (f_oF2) the validity of the International Reference Ionosphere model to predict the time of minimum ionization is checked. Data obtained at different ionospheric stations have been considered. The CCIR and URSI options are used to model calculations. For CCIR option the results show that good predictions were obtained for about 40% of the considered cases. For the rest of the considered data, the model predicts the minimum at times earlier than that observed in the measurements. The percentages of good predictions obtained with URSI option are lower than those corresponding to CCIR one.     

Comparison of ionospheric data from Soviet and Chinese stations with the IRI-86 and the faim models

Median values of ionospheric data from Moscow(55.5° N, 37.3° E), Tashkent (41.3° N, 69.6° E), Beijing (40.0° N, 116.3° E) and Guangzhou (23.1° N, 113.3° E) for different years at LT: 00.00, 06.00, 12.00 and 18.00 were used to produce vertical profiles of electron density with the analytic procedure introduced by Giovanni and Radicella. These profiles were compared with the IRI-86 model profiles adjusted to the F2 experimental peak values and with those obtained with the Fully Analytic Ionospheric Model (FAIM). The main differences with the IRI model are related to the shape of the topside profile under certain conditions and to the presence of the F1 layer in the median values used. The differences with the FAIM model appears to be more marked. The results are discussed.     

全球电离层对2000年4月6-7日磁暴事件的响应

利用分布于全球的电离层台站的测高仪观测数据,对扰动期间,foF2值与其宁静期间参考值进行比较,研究了2000年4月6—7日磁暴期间全球不同区域电离层的响应形态,并通过对比磁扰期间NmF2的变化与由MSISR90经验模式估算的中性大气浓度比(no／nN2)的变化,探讨了本次事件期间的电离层暴扰动机制．结果表明,在磁暴主相和恢复相早期,出现了全球性的电离层F2层负相暴效应．最大负相暴效应出现在磁暴恢复相早期,即电离层暴恢复相开始时间滞后于磁暴恢复相开始时间．在磁暴恢复相后期,一些台站出现正相扰动．研究结果表明,本次事件期间的电离层暴主要是由磁暴活动而诱发的热层暴环流引起的．     

Variability of the ionospheric electron density at fixed heights and validation of IRI-2007 profile’s prediction at Ilorin

A study on the variability of the equatorial ionospheric electron density was carried out at fixed heights below the F2 peak using one month data for each of high and low solar activity periods. The data used for this study were obtained from ionograms recorded at Ilorin, Nigeria, and the study covers height range from 100 km to the peak of the F2 layer for the daytime hours and height range from 200 km to the peak of the F2 layer for the nighttime hours. The results showed that the deviation of the electron density variation from simple Chapman variation begins from an altitude of about 200 km for the two months investigated. Daytime minimum variability of between 2.7% and 9.0% was observed at the height range of about 160 and 200 km during low solar activity (January 2006) and between 3.7% and 7.8% at the height range of 210 and 260 km during high solar activity (January 2002). The nighttime maximum variability was observed at the height range of 210 and 240 km at low solar activity and at the height range of 200 and 240 km at high solar activity. A validation of IRI-2007 model electron density profile’s prediction was also carried out. The results showed that B0 option gives a better prediction around the noontime.     

Validation of B0 and B1 in the IRI 2001 model at low solar activity for Ilorin an equatorial station

A new set of data obtained at low solar activity from Ilorin, Nigeria (geog. latitude 8.5°N, geog longitude, 4.6°E, dip 4.1°S) is used to validate the IRI 2001 model at low solar activity. The results show in general a good agreement between model and observed B0 at night but an over estimation during daytime. The overestimation is greatest during the morning period (0600LT–1000LT). The model prediction for B1 is fairly good at night and during the day. A dependence of B0 on solar zenith angle χ is observed during the daytime. A formulation of the form B0 = A[cos(χ)ⁿ] is therefore proposed. Values of the constants n and A were determined for the period of low solar activity for this station.     

On the use of NeQuick topside option in IRI-2007

The latest version of IRI includes various options for the computation of the topside electron density profile. One of the possible choices is based on NeQuick model. Its inclusion in IRI has been made transferring all the formulations used in NeQuick model. In details, an Epstein layer function is used to describe the electron density profile and the topside shape is controlled by an empirical parameter, connected to the NeQuick F2 bottomside thickness parameter, B2bot. It is computed also in this IRI topside option in order to maintain self-consistency with its original formulation. This paper analyses the possibility of using the IRI bottomside parameters for this option and its impact on the profile and TEC. The case of experimental peak values given as input is also analysed.     

Global validation of IRI TEC for high and medium solar activity conditions

The IRI model offers a choice of options for the computation of the electron density profile and electron content (TEC). Recently new options for the topside electron density profile have been developed, which have a strong impact on TEC. Therefore it is important to test massively the IRI and the new options with experimental data. A large number of permanent stations record dual frequency GPS data from which it is possible to obtain TEC values. Thirty-one worldwide distributed stations have been selected to investigate the capabilities of the IRI to reproduce experimental TEC. Data for years 2000 (high solar activity) and 2004 (medium solar activity) have been analyzed computing modeled values with the IRI-2001 and the IRI-2007-NeQuick topside options. It is found that IRI-2007-NeQuick option generally improves the estimate of the slant TEC, especially in the case of high latitudes stations during high solar activity.     

Comparing TOPEX TEC measurements with IRI predictions

TEC values obtained from TOPEX satellite were compared with the International Reference Ionosphere (IRI) 2001 model estimates. The present work also shows results of the IRI model with the option of a new topside electron density distribution (NeQuick model). TOPEX TEC measurements, which include years of high and middle to low solar activity (2000 and 2004), were analyzed by binning the region covered by the satellite (±66°) every five degrees of modip. In general, there is good agreement between IRI predictions and Topex measurements. Cases with large disagreements are observed at low and high latitudes during high solar activity. Comparing the model predictions using the default IRI2001 model and the NeQuick topside option show that the default IRI 2001 version represents the observed data in a more realistic way, but appears to be less reliable at high and low latitudes in some cases.